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We propose a nanomechanical detection scheme for Majorana bound states, which have been 
predicted to exist at the edges of a one-dimensional topological superconductor, implemented, for 
instance, using a semiconducting wire placed on top of an s-wave superconductor. The detector 
makes use of an oscillating electrode, which can be realized using a doubly clamped metallic beam, 
tunnel coupled to one edge of the topological superconductor. We find that a measurement of the 
nonlinear differential conductance provides the necessary information to uniquely identify Majorana 
bound states. 
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I. INTRODUCTION 

Over the past few years, nanomechanics and the field of 
topological condensed matter systems have been pushing 
limits in their respective area of research. Nanomechan- 
ical systems have proven to be exceptional measurement 
devices for, e.g., mass, force, and positiorPH^, as well as 
a unique platform for studying fundamental questions 
concerning the quantum nature of macroscopic objects, 
theoretically as well as experimentallj^l^. 

Majorana fermions are among the most intriguing fea- 
tures of topological states of matter. They are their own 
anti-particles, i.e., 7 = 7^, and satisfy fermionic anticom- 
mutation relations {7i,7j} — 2Sij. A Majorana fermion 
has half the degrees of freedom of a Dirac fermion. This 
can be seen, for instance, by expressing two Majorana 



fermions ^l,r as 



and 7l — —i{c^ — c). 



where c and are the annihilation and creation oper- 
ators, respectively, of a single Dirac fermion and satisfy 
{c^,c} = 1. Various proposals have been made about 
how to generate Majorana state^i2Hni_ jn view of future 
applications, these states might be particularly useful in 
the field of topological quantum computatiorpS. Unfor- 
tunately, to date no experimental realization of Majo- 
rana fermions has been achieved. Proposed detection 
schemes are based on tunnel setup^S^, interferometer 
setup^^mUl ^j^j ^jjg Josephson effec1p2liniill_ However it is 

fair to say that true qualitative experimental signatures 
of Majorana fermions that persist in realistic systems are 
rather difficult to predict. 

In this paper, we present a novel scheme for detecting 
Majorana bound states (MBS) at the edges of topologi- 
cal superconductors. Our proposal involves an oscillating 
electrode tunnel coupled to a topological superconductor 
(TS) hosting MBS (see Fig.[l]). We show that the inter- 
play between the two weakly coupled MBS, located at 
the edges of the TS, and the oscillating electrode gives 
rise to unique features in the differential conductance of 
the setup. We find that in the presence of the resonator, 
satellite peaks appear in the differential conductance. We 
identify the underlying transport processes giving rise to 



the rich structure in the differential conductance. Fur- 
thermore, we study the dependence of the differential 
conductance on the effective temperature of the oscilla- 
tor and find for an oscillator close to its quantum me- 
chanical ground state a qualitative transport signature 
which is due to the interplay between the resonator and 
the Majorana bound states. 

The paper is organized as follows. We introduce the 
setup under investigation and the underlying model in 
Sec. [nj followed by details on the calculation of the cur- 
rent in the setup in Sec. |III[ In Sec. |IV[ we present our 
results with subsections focusing on a setup with and 
without the resonator. Finally, we conclude in Sec. [Vl 



II. SETUP AND MODEL 

The detection scheme we propose is schematically de- 
picted in Fig. [1] The TS can, for instance, be realized 
as a semiconducting wire with strong spin-orbit coupling 
placed on top of an s-wave superconductor in the pres- 
ence of a magnetic fielcP^t^^. Its left and right edges host 
two MBS which we call 7^ and "fR. One of the edges, say 
the right one, is coupled by tunneling to a movable dou- 
bly clamped beam. While this setup can be technically 
challenging to realize, one should note that the fabri- 
cation of metallic doubly clamped nanoresonators with 
very high freq uenci es and quality factors has been re- 
cently reportecP^I^ni, xhe motion of the beam modulates 
the tunnel amplitude, such that it depends on the dis- 
placement X of the beam. The oscillating electrode is 
held at a bias voltage V and the TS is grounded. 

We approximate the tunnel amplitude as linear in the 
displacement of the beam, ~ ^x^- This approximation 
is justified if the oscillation amplitude is small compared 
to the mean distance between the beam and the edge of 
the TS. The Hamiltonian of the system is then given by 



H — 



H 



MBS 



(1) 



where H^es — J2k^(^)'^k^k describes a spinless electron 
reservoir in the metallic oscillating electrode. Here, we 
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FIG. 1. (Color online) Schematic of the setup. A TS is as- 
sumed to be realized as a ID semiconducting wire on top of 
a grounded s-wave superconductor (SC). The wire can host 
MBS at its left (7^) and right (7^) edges. We assume that one 
of the edges is tunnel coupled to a movable, doubly clamped 
beam (at bias voltage V). The gate electrodes (Vg^-Vg^) can 
be used to increase or decrease the overlap ^ of the MBS by 
changing the effective length L of the TS. 



restrict ourselves to one spin channel only. This is ap- 
propriate since the one-dimensional (ID) semiconducting 
wire on top of an s-wave superconductor in the presence 
of a magnetic field has the same degre es of freedom as 
a spinless Px + ipy superconductop^^'^. In reality, both 
spins will participate and their signal will add up to the 
one we calculate below. Furthermore, we assume a lin- 
ear dispersion e{k). Next, Hqsc = /{2m) + m^l'^x^ /2 is 
the usual harmonic oscillator Hamiltonian describing the 
motion of the beam with an effective mass m and reso- 
nance frequency f2. The term i?MBS = i^lLln/^ char- 
acterizes the overlap between the MBS on the left edge 
and the right edge of the T^^D, Describing these two 
Majorana fermions as one Dirac fermion, we can rewrite 
-f^MBS — ^c^c. The Hamiltonian i?MBS is formally iden- 
tical to that of a single resonant level (RL) at energy 
^. However, due to the nonlocal nature of the MBS, the 
overlap ^ depends exponentially on the effective length L 
of the TS. This distinguishes the MBS Hamiltonian from 
a RL. In order to probe the length dependence of ^ . gate 
electrodes can be installed in proximity to the TS^^ESl^ 
schematically shown in Fig. [T] Then, by applying gate 
voltages, the effective length L, and thus ^, can be tuned. 

The bare tunnel Hamiltonian fftun,o has been intro- 
duced for a related setup in Ref . [11! and the x-dependent 
term i/tun,x is new here. Assuming that the tunneling 
takes place locally at y = 0, both are given by 

Htnnfl = to [^^ iy^O)~^{y^ 0)]jr, (2) 

^ftun,x = -t^£[^Hy - 0) - ^{y - 0)]7fl. (3) 

Note that, due to the form of i?tun,o and i?tun,x, the 
Majorana fermion couples to lead electrons as well as 
lead holes. We cast Eq. ^ into a form containing Dirac 
fermions: 



and analogously for -fftun,x- We assume that the lead 
electrons only couple to jji in -fftun,o and -fftun,x, which 



is justified if the wire is much longer than the super- 
conducting coherence length ^sc of the s-wave SC, the 
characteristic localization length of the MB^^^. 

Since MBS are subgap states, we investigate the regime 
of a large superconducting gap Age, i.e., potQ,£,,eV <C 
Asc, where po is the density of states of the metallic lead. 

If the electrode does not oscillate, the system is de- 
scribed by the quadratic Hamiltonian Hq — Hj-cs + Hosc + 
Hmbs + ^^tun,o- In this case, all Green's functions (GF) 
involving -0 and c operators are known exactly. Since Hq 
does not conserve fermion numbers, the anomalous GF 
do not vanish, e.g., (c(i)c(0))o 7^ 0. Then, nonequilib- 
rium transport properties like the current or the current 
noise can be determined exactljB^. 



III. CURRENT CALCULATION 

We shall treat the i-dependence of the tunneling am- 
plitude as a perturbation to the Hamiltonian Hq. This 
requires that txXzpf\/'2{n) + 1 <^ to, where x^pt = 
l/^/2mU, is the amplitude of the zero point fluctuations 
and (n) the mean phonon number of the resonator. We 
restrict ourselves to second order perturbation theory in 
-fftun.x- Our main focus is the calculation of the cur- 
rent. Putting h — \, the current operator is given by 
/ = -edtN = -ie[H,N], where N = J dy'ip\y)%p{y) de- 
notes the number of fermions in the lead. Using the vec- 
tor notation ^ = (c, c^, V'(0)i '0^(0)) ^ we find I = I0 + /x, 
where 



Here, B is a real 4x4 matrix 



(4) 
(5) 



B 



10 

10 



-1-10 



(6) 



We introduce the fermion GF 



G^^^,{t,t') = G,kM = -i{Tc-^jmk{t'))o (7) 

where j, k £ {1, 2, 3, 4} and Tq denotes the time ordering 
operator on the Keldysh contour. Later we will intro- 
duce Kelysh indices {— ,+} referring to the lower (time- 
ordered) and upper (anti-time-ordered) branch of the 
Keldysh contour, respectively. The average is taken with 
respect to the ground state of the unperturbed Hamil- 
tonian Hq. For ix = 0, the time-independent average 
current can be written as 



do) = 



-f^/£{[G-tM-G-,(..)] 



;^tetH-G^cM]}, (8) 
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where 



(9) 



is the Keldysh Green's function, for which in general we 
have G^{t,t') = G-+{t,t') + G+~ {t,t'). 

Next, we consider the corrections to this current for 
small <x- Using a similar matrix notation allows us to 
write 



with 



/O 1 

10 



V 1 1 0. 



(10) 



(11) 



Introducing the unperturbed oscillator Green's function 

D(i,t') = -*(Tci(t)i(t'))o, (12) 

the first order correction to the fermion Green's function 
can be expressed as (for j, k e {1, 2, 3, 4}) 



/oo 
dri^(-7) 



X J2 A^nD-''{t,r,)GYi{t,r,)Gll{T,X), (13) 

m,n— 1 

where A = A — A''" and a,i3,j E {^■+} denote the 
branches of the Keldysh contour. 

We express the second-order correction to the GF ma- 
trix in terms of the advanced, retarded, and Keldysh 
components as 

G,k{t, t') = U G,k[t, t') Ut = f ^j, % \ (t, i') , 



where the transformation is given by 
1 /I -1 



U = 



V2VI 1 



(14) 



(15) 



This leads to the following compact form of the second- 
order correction to the GF in the rotated Keldysh space: 



*Gf,^(^,t') 



(16) 



dTx dT2 



G(<,ri)AS(Ti,r2)AG(T2,t') 



jk 



where the effects of the oscillator are contained in the 
self-energy 



^fk + ^' ^jk 



(17) 



1 / D^Gt + D^Gf, + D'^Gf^ D"Gf, + D^Gg 



D^Gf,+D^G% 



jk 



We will assume that the mechanical oscillator has a 
very high quality factor, such that the linewidth of the 
uncoupled oscillator is small compared to the effective 
linewidth of the fermionic level at ^. With these assump- 
tions, we can use the following advanced, retarded and 
Keldysh GF in Fourier space: 

D"{uj) = iTx[6{uj + 5{uj - VL)]/[2mU] , (18) 
D^{uj) = -iTi{x^) [5{oj + 5{w + n)] , (19) 



(20) 



with 



p / {m ). Finally, the average current 



including the oscillator can be written as (/) = (/q) -I- 
(-^x,i) + ('^x,2), where (/q) is given by Eq. ([s]) and the two 
remaining terms are (j e {1,2}) 



du) 
2^ 



(21) 



x{[G^)f (0.) - G^f (.)] + IG^'^iu) - Gffi.)]} . 

The analytic result for (/) is too long to be displayed 
here. We demonstrate in the following that the aver- 
age current contains unique information about the MBS 
which can be most easily identified in the nonlinear dif- 
ferential conductance d{I)/dV . 



IV. RESULTS 



Compared to earlier proposals^, we have two addi- 
tional energy scales involved in the detection scheme: 
the resonance frequency 51 of the oscillator and its ef- 
fective temperature Teff. Both are to some extent exper- 
imentally tunable. Assuming that the oscillator is in a 
thermal state, one has {x^)th = 2a;^pf(2(n) -I- 1) where 
(n) = 1/ [exp(r2/roff ) — 1] is the mean phonon number of 
the oscillator and where we set k-Q = \. In the following, 
we will discuss the regime where (n) is small, i.e., compa- 
rable to 1. This is challenging to realize experimentally. 
Note, however, that for the highest resonance frequencies 
of the doubly clamped beams in Ref. [30] (~ 500 MHz), 
the thermal occupation number could actually become 
less than 1 at typical dilution refrigerator temperatures. 
At higher temperatures or lower frequencies, one would 
have to implement additional cooling schemes to bring 
the oscillator to the quantum regime. In our model, we 
make the assumption that the electronic system is at zero 
temperature, T^i = 0. The side peaks in the differential 
conductance can be resolved as long as T^i <C 51. If the 
electronic system and the oscillator are in equilibrium 
with the same heat bath, this translates to the require- 
ment (n) ^ 1. However, T^i <^ 51 can also be fulfilled 
for larger (n) if the oscillator is put into an excited state 
using experimental techniques as, for instance, done in 
Ref.® 
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A. Reviewing the case without the resonator 

For clarity, let us start with the case without oscillator 
(tx = 0) and briefly discuss the dependence of the differ- 
ential conductance on the length of the TS and thus its 
dependence on ^. In the case of a finite overlap of the 
two MBS, the differential conductance shows two peaks 
at eV = ±1^1 which can be seen in Fig. [2] (a) (see also 
Ref. |34| for comparison). Both of the peaks have the 
height 2e^//P^. Importantly, for any finite ^, d{I)/dV at 
V = vanishes in our model. This is in stark contrast 
to the differential conductance in quantum dots coupled 
to one superconducting and one normal leacf^. Hence, 
it is rather straightforward to distinguish this situation 
from the MBS case. As mentioned above, the Hamil- 
tonian Hmbs is formally equivalent to one describing a 
spinless resonant level at energy ^. Therefore, a compar- 
ison to a resonant level case is of pedagogical value. Two 
resonant levels with energies ±|^| might lead to a simi- 
lar signal in the differential conductance as in the MBS 
case. This could, for instance, happen if the magnetic 
field destroyed the superconductivity in an experimental 
realization of our proposal. Increasing the effective length 
of the wire using the gate electrodes, one can tune ^ close 
to zero, which would yield a single peak at V = (not 
shown in Fig.[2|. However, this feature in the differential 
conductance could also be due to a single RL with energy 
^ = and consequently an ordinary bound state might 
mistakenly be identified as a MBS. We conclude that tun- 
neling into resonant levels could be hard to distinguish 
from tunneling into MBS. However, it is fair to say that 
a measurement of d{I)/dV as a function of the variation 
of the length of the TS could yield a strong signature 
of MBS, even in the case of static leads. Furthermore, 
d{I)/dV = 2e^/h on resonance in the MBS case, which 
might not be the case for tunneling into a RL. Interest- 
ingly, we show below that the oscillating electrode ex- 
hibits additional features in the differential conductance, 
allowing for an unambiguous identification of MBS. 



B. Differential conductance including the oscillator 

The differential conductance in the presence of the os- 
cillator and its dependence on the oscillator's tempera- 
ture is shown in Figs.[2][b) and[2|c) for two temperatures 
corresponding to (n) — and 1, respectively. Due to 
the presence of the oscillator, satellite peaks emerge at 
eV = ± In Fig. [S] we depict the energetically 
allowed tunnel processes for (n) = 0, which explain the 
emerging satellite peaks in Fig. [2] (Note that for (n) = 1 
more processes are possible corresponding to the emis- 
sion of a phonon by the oscillator. These are not shown 
in Fig. Isj) Evidently, for the conventional processes in 
Figs. [3[a), [Sjjc), and|3je), the superconducting conden- 
sate does not play any role. Hence, these processes also 
matter for tunneling into a RL where the condensate in 
Fig. |3] would be replaced by a second lead. The pro- 
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FIG. 2. (Color online) Differential conductance for ^/f2 = 0.2 
in the presence of MBS for the case — (without oscillator) 
as well as for the case tx 7^ for different values of (n). The 
satellite peaks have been enlarged for better visibility. 



cesses in Figs, [sjjb), [sjd), andjsf^f), on the other hand, 
rely on the presence of the superconducting condensate. 
Importantly, all processes in Figs.jsj^a) tojslf) contribute 
together to the rich structure in the d{r)/dV shown in 
Fig.[2j Increasing the oscillator's temperature leads to an 
increase in the heights of the satellite peaks and trans- 
forms a dip at eV = ±|f ± fij into a peak. This will 
be further discussed in Fig. [4] As mentioned above, one 
can argue that in a single RL scenario with level energy 
^ = 0, the signal in the differential conductance could 
not be clearly distinguished from the one stemming from 
MBS. Including the oscillator permits us to unambigu- 
ously distinguish the two cases. Figure |4] shows our key 
result. In that figure, we compare the RL case to the 
MBS case, both in the presence of an oscillating tunnel 
contact. For clarity, we focus on a region around the 
resonant frequency f2 of the oscillator and positive bias 
voltage V . Importantly, an oscillator in its ground state 
{{n) — 0) can only absorb a phonon. In the case of the 
single RL, the crucial tunnel process near the resonance 
at eV = Vi (for an oscillator with (n) = 0) is depicted in 
the right inset of Fig. |4] We see that this process only 
sets in for voltages eV > ft (at zero temperature) since 
the state in the right lead has to be occupied in order to 
transfer the energy fl to the oscillator. As a consequence, 
the differential conductance is only positive for eV > 
(shown as a solid black line in Fig. |4]). The situation is 
very different for the MBS case. Here, we have a sec- 
ond crucial tunnel process depicted in the left inset of 
Fig. |4] Then, the dashed-dotted blue line in Fig. [i] illus- 
trates that, for the oscillator in its ground state and MBS 
present, the d{I) /dV is positive for eV < and eV > fl. 



5 



(a) cT-0 V > ^ 



(c) a^c^'ip V>n + { 




e) gtc-i/jt y<-(n-{) 




(b) 




(d) gtcV; l/>n-{ 




--0-- n - ? 




FIG. 3. (Color online) Schematic of energetically allowed 
tunnel processes for an oscillator in the ground state. Panels 
a)-f) depict the processes stemming from the tunnel terms 
c'lp, c'-ip^ a^c^ip, a'cip, a^CTp^, and a^c'ip^ respectively, 
where x — + a. Black dots (circles) represent electrons 
(holes) and the condensate is depicted as the blue bubble. 



This feature is only due to tunneling through MBS and 
clearly separates the RL from the MBS scenario. If we 
gradually increase (n) from zero to one (where only the 
two extremes are shown in Fig. |4]), the dip for (n) = 
in the MBS case transforms smoothly into a peak due to 
additional processes where the oscillator emits a phonon. 
The dip at eV = is due to an interference effect be- 
tween the two participating tunnel processes which can 
be understood on the basis of Fermi's golden rule. To 
illustrate this, we choose ^ = and (n) — and write 
the system Hamiltonian as H — Hq + i/tun with 



(22) 
(23) 



In our initial state, the Fermi sea should be filled up 
to the chemical potential eV (measured from ^ = 0), 
furthermore the fermionic subgap state c is empty and 
no phonons are present, written as 

N) = |0^,0„0,) , (24) 

with energy Ei = Ep, Ep being the energy of the filled 
Fermi sea. 

To the lowest order in ix, the leading contribution to 
d{I)/dV at eV « $7 is generated by a combination of the 
two tunneling processes c^ip and ^pc, where one of them 
is accompanied by the creation of a phonon mode. 

The current-carrying transport processes we focus on 
are mediated by the processes described by c''ip in combi- 
nation with (cf. insets in Fig. |4]) . For these transport 
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FIG. 4. (Color online) Differential conductance for ^ ~ at 
eV ~ Q for a single RL coupled to an oscillator with (n) = 
(solid black line) and (n) — 1 (dotted gray line). In the former 
case, the oscillator can enhance transport only for eV > SI. 
This is different in the presence of MBS. Then, an oscillator 
with (n) — gives rise to a positive d{I) /dV for eV > fl and 
eV < fl (dashed-dotted blue line). We also show in dashed 
red the d{I) jdV for (n) = 1 in the MBS situation. The insets 
depict the decisive tunnel processes for the (ji) = case. 



processes the final state lacks two electrons in the Fermi 
sea (at momenta |fc|, \'k'\ < kp\ contains one phonon and 
an unoccupied subgap state c, i.e. we have 



with energy 



I/) = V'feV'fc' |0v.,0c, la) 



Ef = Ep- e(k) - e{k') + , 



(25) 



(26) 



where e{k) is the single-particle electron energy and the 
phonon energy is fl. As usual, the transition amplitude 
between the states above, up to second order in fftun is 
given by 



A 



1 



Ei — Hq 



(27) 



and according to Fermi's golden rule, the rate F is given 
by summing over all final states respecting energy con- 
servation 



r = 2nJ2\AM'S{Ef-E,[ 



(28) 



k,k' 



The rate F is proportional to the current, and therefore 
dr/dV proportional to the differential conductance. By 
evaluating the amplitude Afi and the rate F, we find 
that for the above process, the amplitude Afi vanishes 
for tunneling particles with energy close to eV ~ fl. This 
leads to the conclusion that for eV « the differential 
conductance vanishes due to an interference effect be- 
tween the two processes (of order io^x) generated by the 
terms c'ipito + t^a') and 'ipc{to + t^a') in the tunneling 
Hamiltonian. This explains the dip in the differential 
conductance d{I) /dV in Fig. |4] 
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V. CONCLUSION 

To conclude, we have presented an idea of coupling Ma- 
jorana bound states to a sensitive nanoelectromechanical 
measurement device. We have shown that a setup where 
an oscillating, doubly clamped beam is tunnel coupled to 
a topological superconductor gives rise to unique trans- 
port signatures based on the interplay between the me- 
chanical excitations and the Majorana bound states. A 
smoking gun signature of Majorana bound states has 
been identified for an oscillator close to the quantum 
ground state which can be achieved by cooling a 500 MHz 
resonator to T 20 mK. This energy scale is well below 
the large parameter of our model, i.e. the superconduct- 
ing gap Asc- When we, for instance, take an InAs wire 
proximity coupled to an Nb s-wave superconductor, the 



superconducting gap of Nb AgQ ~ 15 K can induce a su- 
perconducting gap of Ag^A** ~ 1 K in the wir^^^. Hence, 
our predictions are in principle observable at dilution re- 
frigerator temperatures. 
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